Exact solution of the infinite-^/ Hubbard problem and other models in one dimension 
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An interacting spin-fermion model is exactly solved on an open chain. In a certain representation, 
it is the nearest-neighbor Hubbard model in the limit of infinite U (local interaction) . Exact solution 
of its complete energy eigen-spectrum is accomplished by introducing a unitary transformation which 
maps the original problem to a tight-binding model of the fermions only. Physically, the exact 
solution implies the absence of Nagaoka ferromagnetism in the ground state for arbitrary electron 
densities. The present method solves a class of very general models exactly. Few more problems are 
discussed as an application of this unitary transform method. 

PACS numbers: 71.10.Fd,71.10.Hf,75.10.Lp,05.30.Fk 



I. INTRODUCTION 

The physics of strongly correlated lattice electrons is 
complex as well as interesting. The electronic properties 
of very many exciting materials (such as, the transition 
metal oxides, rare earths etc.) exhibit this physics in a va- 
riety of different ways. A minimal model for understand- 
ing these systems requires one to consider, at least, the 
local coulomb repulsion, in addition to the tight-binding 
hopping of electrons. In an effective one-band system, 
the local repulsion, U, is the energy cost for putting two 
electrons with opposite spins together on the same site 
(Wannicr orbital). This model, famously called the Hub- 
bard model, has been a subject of great interest. It poses 
one of the most difficult problems in quantum many body 
theory. Although the model has been long solved ex- 
actly in one dimension (ID) using Bethe ansatz 1 , and 
also shown to be integrable 2 , it still evades an exact an- 
alytical solution in higher dimensions. In recent times, 
the interest in the Hubbard model has been further re- 
inforced by the high-Tc superconductivity in cuprates, 
and also by the need to understand other phenomena in 
strongly correlated electrons 3 . 

There has been a long-standing interest in under- 
standing the metallic-ferromagnetism through Hubbard 
model. The limit of infinite local repulsion in one-band 
Hubbard model presents an interesting case study in this 
context. In this limit, the ground state on certain lat- 
tices (for example, square lattice) was shown to be satu- 
rated metallic-ferromagnetic for a single hole in an oth- 
erwise half- filled system (Nagaoka- Thouless theorem) 4,5 . 
Subsequent variational and numerical studies have shown 
that the Nagaoka ferromagnetism survives for finite hole 
densities away from half-filling, up to a (lattice depen- 
dent) critical doping 6 ' 7 ' 8,9 . The Lieb-Mattis theorem, 
however, rules out the existence of ferromagnetism in 
ID 10 . For the infinite-C/ Hubbard model with nearest- 
neighbor hopping, it implies the lack of Nagaoka ferro- 
magnetism, which is borne out by the analytic studies 
of this problem for one and two holes explicitly. These 
investigations either use the Bethe ansatz approach 11 , or 
work with an effective spin Hamiltonian in the presence 
of a single hole 12 , starting with a Gutzwiller projected 



Hamiltonian. 

Recently, we have developed a new approach to the 
infinite-?/ Hubbard problem 13 . In our formulation, wc 
canonically represent an electron in terms of a spinlcss 
fermion and the spin-1/2 (Pauli) operators. We then 
write the Hubbard model in this representation. Finally, 
we take the limit U — ► oo, and get the infinite-?/ Hamilto- 
nian, . Although our prescription is applicable to any 
lattice, it clearly distinguishes the bipartite lattices from 
others. By exploiting the two sublattice structure of a 
bipartite lattice, we can write Hoo in a beautiful form 
(resembling the Anderson-Hasegawa Hamiltonian, but 
different from it and fully quantum mechanical) which 
reveals the phenomenon of metallic-ferromagnetism in 
the infinite-?/ Hubbard model in a very transparent way. 
(Otherwise, we all know that the Nagaoka ferromag- 
netism is a non-obvious strong correlation effect.) 

In the present study, we investigate ID infinite-?/ Hub- 
bard problem within our approach, and exactly solve it 
by means of the suitably constructed unitary transfor- 
mation. From this exact solution, wc conclude that its 
ground state is 'correlated metallic' and 'ideal paramag- 
netic' for arbitrary density of electrons. Ours is a 'non- 
Bethe' method which also solves a class of very general 
models in ID. The infinite-?/ Hubbard model happens to 
be one among them. This paper is organized as follows. 
First we present the exact solution of Hryo in ID. Then, 
we identify a general class of models which can be ex- 
actly solved by our method. We also briefly discuss the 
Andcrson-Hascgawa problem and the minimal coupling 
lattice Hamiltonian in the light of ideas developed here. 
Finally, we conclude with a summary. 



II. 



INFINITE-f/ HUBBARD MODEL IN ID 



According to a recently developed canonical (and in- 
vcrtible) representation, an electron can be described in 
terms of a spinless fermion and the Pauli operators 13 . On 
a bipartite lattice, we can represent the electrons on dif- 
ferent sublattices in two different, but equivalent, ways. 
(In principle, we can generate infinitely many equivalent 
representations through unitary transformations, but our 
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purpose is served by the following two forms.) The elec- 
tronic operators in this representation arc: fj^ = ; cr ; + , 
fl i = (i'ipt — 0;crf)/2 on the odd-numbered sites (I = 
1,3,5,...), and = ifa<r+, = - ijn<rf)/2 
on the even-numberd sites, where <fii = (aj + a{) and 
vfe = (a, — &t) are the Majorana fermions. Here, a l and 
a[ arc the spinless fermion operators, and <?/ the Pauli 
operators, on the I th site. Moreover, the electronic num- 
ber operator is written as: N e = L + 52j = i.(l — ^z)°f> 
and the 'physical spin' of an electron on the I th site is 
given by: Si = hiai/2, where hi = aja^ For complete- 
ness, also note the local mapping: |0) = |o), |T) = \t), 
||) = 7), and |t|) = \t) (on an odd-numbered site) and 
= — |o) (on an even-numbered site). Here, |0) , ||), \l) 
and |t|) are the local electronic states (in the usual nota- 
tion). The states |o) and |») denote the empty and filled 
site, respectively, for the spinless fermion, and {|+}, |— )} 
is the basis set of a single Pauli spin. Note that |+) and 
|— ) represent the actual electronic spin on a site when it 
is occupied by a spinless fermion. 

Now, consider the Hubbard model with nearest- 
neighbor hopping on an open chain. (The reason for 
working with an 'open' chain will become clear shortly.) 
Since the hopping process is bipartite, we can use the 
above two forms of the representation for electrons to 
convert the Hubbard model into a corresponding 'spin- 
fcrmion' model. In the limit of infinite-?/, we get the 
following Hamiltonian. 

L-l 

Hoo = ~t^2 X l,l+i (X^i+i + aj+i^i) (!) 
1=1 

For the details of its derivation, please refer to Ref. [11]. 
Here, X l l+1 = (1 + <?; ■ 07+1 )/2 is the Dirac-Heisenberg 
exchange operator, and L is the total number of lattice 
sites. 



A. Exact analytic solution 

For the moment, we forget the physical origin and pur- 
pose of -Hoc 7 and just take it as a given spin-fermion 
model in one dimension. Our immediate goal is then 
to find its eigenvalues and eigenstates. There are two 
key features that we make use of in exactly solving this 
problem. First is the property, Xf l+l = 1, of the ex- 
change operators. And, the second is the open boundary 
condition of the ID lattice (similar to the XY spin-1/2 
chain 14 ). We exploit the former to construct a unitary 
operator which, on an open chain, transforms Eq. (1) to 
a tight-binding model of the spinless fermions only. 

We first define the following unitary operator on the 
first bond [that is, for the pair of sites (1,2)]. 

U1.2 = (1 - h 2 ) + h 2 X 2tl (2) 
Here, X 21 = X 12 . Clearly, li\ 2 = U 12 and U\ 2 = 1. 
Thus, Ui^2 is both Hermitian as well as unitary. More- 
over, it has the following important property. 

li\ 2 (a\x 12 a 2 ^ U 12 = a\a 2 (3) 

In the above equation, leaves a\ and Xi,2 unaffected 
while transforming a 2 — ► A2.1a.2- Thus, Ui >2 gets rid 
of the exchange operator Ai.2, and what remains is the 
hopping of the fermions alone. As it happens, we will get 
rid of the exchange operators on each bond by carefully 
following this approach. 

Before constructing a similar U 2 ^, for the next bond, it 
is important to consider the effect of U\^ 2 on other terms 
in -Hoc,. Clearly, Ui 2 leaves the operators on other bonds 
unaffected, except for the bond (2,3). The net effect of 
this unitary transformation on is the following. 



U\ 2 i/oo ^1,2 = — t I f a\'a 2 + a 2 a 1 ) + (a 2 X\. 2 X 2 ^ a 3 + a 3 X3, 2^2,1 a 2 



it; 



(4 



L-l 

£ 

1=3 



&i -^1,1+1 a i+i 



a 



X 



1+1 ^1+1,1 a i 



(4) 



The Ui t2 transfers Xl,2 from the first bond to the second. 
Wc define ,Y 2 .3 = .Y1.2.Y2.3 and X 3t2 = A : >. 2 A 2 .i. which 

7^ ^2,3 

X 3 2 X 2 3 = 1. Thus, ^2,3 is unitary but 



replace A2.3 on bond (2,3). Note that X 2 3 



and Xl 3 X 23 



not Hermitian (unlike A2.3). 



Now wc define the following unitary operator for get- 
ting rid of X2.3 and X 3 2 from the terms inside the second 



parentheses in Eq. (4). 

^2,3 = (1 - n 3 ) + h 3 X 3 , 2 (5) 
Note that U.2,z is unitary but not Hermitian (unlike U\a)- 
We can show that ^2,3 is invariant under k/2,3, while 



Therefore, 



U 2 3 a 3 U 2 ^ — a 3 X 3 2 . 



^2,3 ("2 ^2,3 "3 ) ^2,3 — 4«3- 



(6) 
(7) 



3 



Moreover, U\ 3 (a\X^,^a A )U 2 3 



0,3^3 40,4, where X% t \ = 

■^2,3^3,4 = ^1,2^2,3^3,4- It is clear by now that we 
can continue this process, and get rid of all the exchange 
operators in H^. To achieve this, we define the following 
general unitary operators, 

= (1 - ni+i) + fn +1 X l+1 j (8) 

L-1 

U = Y[Ui,i +1 (9) 

1=1 

where <£(,;+i = ^-i^-X^+i = Y\l n= iX m . m+1 , for I = 
2,L — 1. We can show that the full unitary operator 
U, of Eq. (9), transforms Hoo to a Hamiltonian of free 
spinless fermions. That is, 

L-1 

U ] H x U = -tJ2 («l«;+i + 4f 1^) (10) 
;=i 

= 5Z £fc & k & k ( n ) 
fc 

where = — 2icosfc, and k is the momentum. The 
operators {a^} arc the fcrmions in the momentum space. 
Since Eq. (10) is derived on a chain with open boundary 
condition, the Fourier transformation between {ai} and 



{a,k} is defined as: &i = y XTT Sfc °* sin/cZ, where k = 
nir/(L + 1) for n = 1, 2, • • • , L. Hence, the exact solution 
of Hoc in ID. To this end, we also note that 



NU — N and M a U = M a 



(12) 



where TV = X)/=i ^ s the number operator of the spinless 
fermions, and M a = Yli=i G i ^ s the total a z operator. 
Therefore, U diagonalizes a more general Hamiltonian: 
H = — ATV — ijM a , where A is the chemical potential 
of the spinless fermions and 77 is the external 'magnetic' 
field acting on the spins. 

The above exercise presents a rigorous and transparent 
case of the complete decoupling of the Pauli and Fermi at- 
tributes of the electron. As a result of this decoupling, the 
energy eigenvalues of Hoo become independent of spins, 
giving rise to an extensive entropy in every cigenstatc. 
Just as the spinless fcrmions {a;} transform under U, 
the spins {(?;} also transform to new spins. However, 
the total a z is invariant under U, as noted in Eq. (12). 
It is these transformed spins that are absent in the free 
spinless fermion Hamiltonian [Eq. (10)]. 



B. Absence of the Nagaoka ferromagnetism in ID 

In order to discuss the ground state of the infinite- 
ly Hubbard model, first let us recall the complete con- 
nection between the problem worked out in the previous 
subsection and the physical Hubbard model. Technically, 
the infinite-?/ Hubbard model in our representation is not 
just -ffoo, but Hoc + 00 — ^0- That is, the physical 



problem is described by H&Q with an infinite chemical po- 
tential for the spinless fermions (please look into Ref. [11] 
for details). Next, we note that TV = N e — 2Nd, where 
N is the total number of spinless fermions (that is, the 
number of singly occupied sites in terms of the electrons), 
N e is the total number of electrons, and is the total 
number of doubly occupied sites. Clearly, TV is a con- 
served quantity of Hqo [evident from Eq. (1)]. So is N e , 
and hence Njj. It is evident from the following explicit 
form of iJoo in terms of the electron operators. 



L-1 



1=1 s 



n l+l,s) 



ft,sfl+l, 



h.c 



(13) 

Since N e = N + 2Nd, we can label the sectors of states 
for a given N e in terms of the partitions: (N, Njj). The 
physical validity of a partition, however, depends upon 
whether N e < L (less than or equal to half-filling) or 
N e > L (more than half- filled case), subject to the nat- 
ural constraint: N + No < L. That is, an arbitrary 
partition of the integer N e into two other integers N and 
Nd does not necessarily denote a physical sector. For 
< N e < L, the constraint is guaranteed to be satis- 
fied. Therefore, all partitions are valid physical sectors. 
For example, if N e = 7 and < L, then the correspond- 
ing sectors of states are: (7,0), (5,1), (3,2), and (1,3). 
However, when N e > L, the constraint will disqualify 
many partitions. For example, if L = 4 and N e = 7, 
then the only physical sector is: (1, 3), as the other par- 
titions such as (3, 2) don't respect the constraint. In gen- 
eral, the physical sectors for N e < L are given by the 
set: {(N e - 2N D , N D ), V N D = 0, 1, • • • , [N e /2]}, where 
[N e /2\ = N e /2 for even values of N e and = (N e - l)/2 
for odd values. For N e > L, we use the relations: N[2L — 
N e ] = N[N e ], andN D [2L-N e ) = L-N e +N D [N e ], to find 
the physical sectors. These relations are a consequence 
of the particle-hole transformation on the electronic op- 
erators. Here, N[N e ] and Nu[N e ] denote the dependence 
of N and N D on N e . 

The exact solution of Hqo gives a highly disordered 
ground state in terms of the spinless fermions and the 
Pauli spins. However, we need to carefully translate its 
meaning for the electrons. Interestingly, we are able to 
show that the ground state of the infinite-C/ Hubbard 
model for a given N e , is a fermi-sea which is 2 N -fold 
degenerate [where N* is defined in Eq. (14)]. Physically, 
it means that the ground state is metallic and ideally 
paramagnetic. In other words, it is not Nagaoka ferro- 
magnetic. Nor it is a kinematic singlet (like a normal 
electronic fermi-sea). 

Now, the proof. Due to the fact that U = 00 and 
it is the chemical potential of the spinless fermions, the 
ground state of the Hubbard problem, for a given 7V e , 
lies in the sector (N*,(N e - N*)/2), where TV* is the 
maximum allowed value of TV for the given TV e . 



TV* = 



TV,. 
2L 



N, 



< TV e < L 
L< N P <2L 



(14) 
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Since N is invariant under U [Eq. (12)], the infinite-?/ 
ground state corresponds to the fermi-sea of TV* spinlcss 
fermions with dispersion ek [Eq. (11)]. We derive the 
following exact expression for the ground state energy. 



TT N* 

Eg [N*] = -2tcos[- — 



sin 



TT N* 

2 L + l 



2 L + l 



(15) 



Expectcdly. E g — for N e = and 2L (trivial cases: 
empty and fully filled bands, respectively), and also for 
N e = L (the half-filled case for U = oo). In the ther- 
modynamic limit (L — > oo). for a finite electron density, 
n e = N e /L, the ground state energy density, e g = E g /L, 
can be written as 15 : 



2t 

e g [n e ] = 1 sin(7m e )|. 

TT 



(16) 



Now, we enumerate the spin degeneracy in the ground 
state, which will decide the magnetic nature of the 
ground state. Let us discuss N e < L case first. In this 
case, the ground sector is (N e , 0). That is, the number of 
electrons is completely exhausted by the number of spin- 
less fermions. As No = 0, the remaining L — N e sites 
must be empty. The state of an empty site is uniquely 
|o). However, the state of a site, occupied by a sin- 
gle electron, is either \t) or |^), corresponding to the 
fact that it could be an f or J. spin electron. Therefore, 
corresponding to any given distribution of N e spinless 
fermions on L sites (any one of the l Cn c combinations), 
there are exactly 2 N * states (a total of l Cm c x 2 n " states 
in the ground sector). For example, on a chain with 
L = 7 and N e = 3, the states in the ground sector are 
like: jJJJoooo), where the filled sites could be in any 
one of the 7 Cs combinations. On each filled site, the 
Pauli spin could be + or — (without affecting the num- 
ber of electrons). Hence, 2 3 different M a states. Coming 
back to the general situation, these 2 Nc states can be 
grouped according to their M a values { N "CM a states for 
a given M a ). We know that M a is a conserved quan- 
tity of Hoo, and it is also invariant under U [Eq. (12)]. 
Therefore, 2 N " different M a states in the ground sector 
will be degenerate, as the exact energy eigenvalues of 
Hoo are independent of M a . Hence, the ground state, 
in the sector (N e ,0), is a 2 JVe -fold degenerate fermi-sea 
of N e spinless fermions. Since there are only empty or 
singly occupied sites, the 2 JVe -fold degeneracy is strictly 
due to the physical spin of electrons. Therefore, we con- 
clude that the exact ground state of the infinite-t/ Hub- 
bard model is ideally paramagnetic and metallic (more 
correctly, strange or correlated metallic, as it is not a 
fermi-sea of the normal electrons). 

Furthermore, we note that the arguments for N e > L 
are the same as that for N e < L. The (only) key dif- 
ference between the two is that the for N e > L, there 
are only singly or doubly occupied sites in the ground 
sector, while for N e < L the sites are either empty or 
singly occupied. For example, when N e = 11 and L = 7, 



then a typical state will be of the form: | . . . ft 1 1) ■ 
This state is the counterpart of a previously mentioned 
state for less than half-filling. The physics of the Hub- 
bard model on a bipartite lattice for N e electrons is same 
as that for 2L — N e electrons. Without going into the 
(repetitive) details of the analysis all over again, we con- 
clude that there is no Nagaoka ferromagnetism in the 
one-dimensional inifinte-[/ Hubbard model, and that the 
exact ground state is a 2 N -fold degenerate fermi-sea of 
N* spinlcss fermions, where N* is defined by Eq. (14). 

To this end, we would like to make four comments. 
First is that the state labeling and counting procedure 
presented above is applicable to all lattices. It is not spe- 
cific to ID (although the solution is). Second comment is 
that, in the cases different from the present ID problem, 
it will be impossible to completely get rid of the exchange 
operators. Due to which, the different M„ states are not 
guaranteed to be degenerate. Therefore, we stand a clear 
chance of finding some sort of metallic magnetism (ferro 
or antiferro 12 or something else) on other lattices (and 
hopping geometries). Third is a minor comment about 
U = — oo problem. In this case, the ground sector corre- 
sponds to N D = N e /2 (for even N e ) or (N e - l)/2 (for 
odd N e ). The ground state (say for even N e ) is a L Cm d - 
fold degenerate hard-core bosonic state with no kinetic 
energy gains (due to spinless fcrmion). The final com- 
ment is about the finite temperature calculations. Since 
U is infinite, the ground sector (for a given N e ) is the 
only part of the Hilbcrt space which is accessible by finite 
temperatures. The thermodynamics of this problem can 
therefore be worked out in the canonical ensemble of the 
spinlcss fermions. The Pauli spins remain ideally para- 
magnetic down to absolute zero temperature. This sets 
the quantum coherence temperature for electrons to be 
zero (even though there is this Fermi temperature scale 
for the spinless fermions) 16 . 



C. Spin-orbital model 

In a different incarnation, the Hamiltonian can 
be considered as a one-dimensional model of the coupled 
spin and orbital degrees of freedom. By applying the 
Jordan- Wigncr transformation on the spinless fermions, 
we can derive the following Hamiltonian, 



L-l 



H, 



so 



'+i vi T i+i i ■( i+i/ 9 
1=1 " I 

a (17) 

where {a/} have been changed to {r ; }, and r ; z = 2fii — 1. 
In a transition metal ion with two-fold orbital degener- 
acy, the orbital degrees of freedom can be described in 
terms of the Pauli operators. Let f denote this orbital 
degree in the present discussion. The model Hso is spe- 
cial case of the more general Kugel-Khomskii type mod- 
els 17 . Equation (17) thus presents an exactly solvable 
spin-orbital lattice model, in which the spins and the or- 
bitals behave as decoupled. While the orbital part acts 
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as an XY chain, the spins become paramagnetic. 

III. GENERAL CLASS OF EXACTLY 
SOLVABLE MODELS 

While transforming to the tight-binding model of 
spinlcss fcrmions, it became clear that a very general 
class of models can be solved exactly by our method. A 
Hamiltonian in this class can be written as: 

L-l 

H = -t E + «!+i^+i.i«z) ( 18 ) 

1=1 

Here, Ti,i+i is some unitary operator on the bond (7, l+l), 
and by defintion Ti+ij = T^ l+l . Moreover, T;,z+i doesn't 
have to commute with Tj-i^ and while strictly 

commuting with the T operators on other bonds, and 
also with the fermions, {a;}. With a few careful steps of 
algebra, this Hamiltonian can be transformed to a tight- 
binding model of the fermions [Eq. (10)] with the help of 

U = rL=i Ul,l+U where Uij+i = (1 - fn +1 ) +ni + iTi + ij. 
Here, = n^i^.m+i and li+ij = For 

a special case in which T;,z+i are just the phase factors 
e KM+i ; we can g e t r j(j f these on a Caley tree of arbitrary 
coordination, z (for a nearest-neighbor chain, z = 2). In 
this case, it can be done not only for the fermions but 
also for the bosons. 



for higher spins. Since Xf 2 = 1, just like in the spin-1/2 
problem, the corresponding spin-fcrmion model [that is, 
Eq. (1)] can be diagonalized in the same way. 

B. Some physical corrolaries 

1. Anderson-Hasegawa problem 

Our method of getting rid of the unitary fac- 
tors has an interesting consequence for the Anderson- 
Hasegawa (AH) problem. The AH model, Hah = 

-*E(t.m> v" 1+ °2 ( & W* l,m& rn. + h.c), describes the 
motion of locally spin-projected electrons on a lattice, 
with classical spins, {$"2;}, in the background. Very of- 
ten, the phases <£>; !m arising due to the spin-projection 
along fli and n m are ignored without proper justifi- 
cation. There have been studies which rightly empha- 
size on taking into consideration the effects of these 
phases while computing physical properties 18 . The spe- 
cial case on a Caley tree discussed earlier, however, im- 
plies that indeed the original Haf can be transformed 

to H AF = -t£<j,m) \/ 1+C1, 2 nm (ata m + h.c.) on Caley 
trees. Thus, we give a reason for dropping the phases in 
Haf, at least on a Caley tree. On an arbitrary lattice, 
however, one must keep these phases. 

2. Minimal coupling Hamiltonian in ID 



A. Spin-fermion model for higher spins 

As an academic exercise, we construct the spin-fermion 
models for higher spins, which belong to this general class 
of exactly solvable models. We achieve this by construct- 
ing the 'correct' analog of the exchange operator for a 
pair of higher spins. By correct exchange we mean that 
l^ii)!? 7 ^} mus t become \m2)\mi) under the exchange op- 
erator. For clarity, we work it out explicitly for the spin-1 
case. Here, S z \m) — m\m), with m — 1,0,1 (-1 is de- 
noted as 1; for spin-S 1 , m = S, S + 1, . . . , S). The usual 
spin-spin interaction, Si • S2, does not really exchange 
mi and 7712. Hence, we construct an operator such 
that Xi i2 \mi}\m 2 ) = |m2)|mi). For example, |0)|0) re- 
mains unaffected under Xi j2 , while |1)|1) becomes |1)|1) 
and vice versa. Explicitly, in terms of the spin operators, 
this spin-1 exchange operator can be written as: 



X 



1 



1 



1,2 



= 1 - (St, + S 2 2z ) + -Sx • S 2 + - (S x ■ S 2 ) 



(Si • S 2 )(S lz S 2z ) + -(Sj x S 2 ) 2 (5 lz - S 2z ) 



-2KS1 xs 2 ) 2 ] 2 



(19) 



where (Si x S 2 ) z = ^(Si + S 2 - - Si-S 2+ ). It is clear 
that we can similarly construct the exchange operators 



We now briefly discuss the Peicrls minimal coupling 
of the quantized electromagnetic radiation to the lattice 
fermions, in the light of these gauge removing tricks. 
Consider tight-binding electrons with nearest neighbor 
hopping on an open ID lattice. The corresponding gauge- 
invariant Hamiltonian can be written as: 

H = -* E E (&i,X e< * ir A * dx + h - c ) + H ? 

1=1 s 

(20) 

Here, e = — |e| is the electronic charge, Hp = 
E q a ^(^qA^qA + \) ls the Hamiltonian (where 

a q A, a qA are the Bose operators for an electromagnetic 
field of wavevector q and the polarization A), and A x is 
the x-component of the vector potential. We have cho- 
sen the x-axis to be along the chain (y = z = line). In 
general, the vector potential can be written as: 



AM=iYj— — 



"qA(r)d qA 



«qA( r )a q A (21) 



where u qA is a normal mode (vector) function (including 
the information on polarization). Clearly, the vector po- 
tential operators at different spatial points commute with 
each other. The corresponding electric field operator is 
given by: E(v) = i[l(r), Hp]. 
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Now again, we get rid of the field dependent factors 
from the hopping by applying the following unitary trans- 
formation: 



U = cxp < i 



i^r^hl / A x dx > 
1=2 ) 



(22) 



where hi = hi^ -\-hi,i- Under this U, we get the following 
transformed Hamiltonian: 

L-i hi 

U^HU = -tJ2 E(#+i,X + h-c)+U^H F U (23) 

1=1 s 

While the hopping becomes simple, the field Hamiltonian 
transforms to: 

U^H F U = H F -eY]hi / E x dx + 

1=2 Jl 



' l,l'=2 



(24) 



Here, the second term is the potential energy in the 
presence of electromagnetic field. Suppose E x is the 
electric field of a free radiation propagating in the z 
direction. Then, e^^ 2 n; J 1 E x dx is same as P X E X , 

the dipole interaction. Here, P x = e^2 l=2 (l ~ i s 
the electric polarization operator. Furthermore, V/.z< = 

y Ji dx Ji dx' X) q A ^-i u x qX u x'.q.x] i s the 'Coulomb' repul- 
sion between electrons, generated by the 'exchange' of 
the photon (facilitated by U). Thus, we have derived the 



gauge independent Coulomb and the dipolc-ficld inter- 
actions, starting from the minimal coupling Hamiltonian 
on a lattice, without any approximations. Since the vec- 
tor potential commutes at different points, in principle, 
we can do the same on a Caley tree as well. 



IV. 



CONCLUSION 



To summarize, we have exactly solved the infinite-C/ 
Hubbard model with nearest neighbor hopping on an 
open chain. We use a newly developed canonical rep- 
resentation for electrons, in which the Hubbard model 
becomes a spin-fermion model. This spin-fermion model 
is exactly solved by applying a non-local unitary transfor- 
mation. Under this transformation, the Pauli spins com- 
pletely decouple from the fermions, as a result of which, 
the ground state is correlated metallic and ideal param- 
agnetic for arbitrary density of electrons. 

This method solves a class of very general models. 
Guided by this observation, we have also constructed 
spin-fermion models for higher spins, by suitably extend- 
ing the notion of 'exchange' operators for higher spins. 
It is explicitly worked out for spin-1. (The spin only 
models, using our definition of the exchange operator 
for higher spins, exhibit interesting properties. These 
calculations will be discussed elsewhere.) By using the 
ideas developed here, we have shown that the phase fac- 
tors in the Anderson-Hasegawa model can be droped on 
Caley trees. We have also derived the 'dipole' interac- 
tion and Coulomb repulsion starting with Peierls minimal 
coupling Hamiltonian. 
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